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Abstract 

Long ago appeared a discussion in quantum mechanics of the prob- 
lem of opening a completely absorbing shutter on which were imping- 
ing a stream of particles of definite velocity. The solution of the prob- 
lem was obtained in a form entirely analogous to the optical one of 
diffraction by a straight edge. The argument of the Fresnel integrals 
was though time dependent and thus the first part in the title of this 
article. In section 1 we briefly review the original formulation of the 
problem of diffraction in time. In section 2 and 3 we reformulate 
respectively this problem in Wigner distributions and tomographical 
probabilities. In the former case the probability in phase space is very 
simple but, as it takes positive and negative values, the interpretation 
is ambiguous, but it gives a classical limit that agrees entirely with our 
intuition. In the latter case we can start with our initial conditions in 
a given reference frame but obtain our final solution in an arbitrary 
frame of reference. 



1 Introduction 

Long ago |J one of us (M.M.) discussed in quantum mechanics the problem 
of opening at time t = a completely absorbing shuttler situated at x = 0, 
on which were impinging a stream of particles of definite velocity. In units 
in which h and the mass m of the particles are unity, the problem reduces to 
finding a wave function that satisfies the free one dimensional time dependent 
Schroedinger equation i.e., 



.dijj(x,t) 



ld 2 ?p(x,t) 
2 



dt 2 dx 2 

with the initial condition 

ip(x, 0) = exp(ikx)0(— x 

where 9(x) is the step function 

e( x ) = { 1 if x>0 

1 ' to if x<0 



(i) 



(2) 



(3) 



The solution of this problem was given in reference 1 and later Nussensveig 
gave it the name M(x, k,t) and it can be expressed as[l,2,3] 



M(x, k,t) — - exp 



i(kx — -k 2 t 



erfc(e~ l77/i w) 



e m / 4 exp 



i{kx k 2 t) 



1 

V2 



- 2 -c( w ) 



+ i 



2 - S{m ' 



where 



x — kt 



w 



2t 



and the error integral is 



2 f 00 2 

erfc(z) = —= \ e~ y dy, 

\/TT Jz 



(4) 



(5) 



(6) 



while the Fresnel integrals are defined by 



C(w) 



IT JO 



cosy 2 dy, S(w) 



TT JO 



siny 2 dy. 



(7) 



Although we have assumed k real, as in the units we use it is the velocity 
or momentum of the impinging particles, all the above expressions remains 
valid for complex k so long as Im k < 0. In that case we have the alternative 
representation p], |3j 

M(x,k,t) = — -^ dre, (8) 

which follows from the fact that both sides are solutions of ([[]) satisfying the 
initial condition (fj). The Green function for the one dimensional free particle 
Schroedinger equation has the form 

U(x-af,t) = eXp[i{x ^ 2/2t \ (9) 

as it satisfies Eq. (ffl) for any t > 0, but when t = it becomes the delta 
function 5(x' — x). As our initial condition is (0) it is clear[l,2] that our 
function M(x, k, t) can also be written as 

M(x,t) = U(x-x',t)exp(ikx')dx'. (10) 

J — oo 

The expression 

\M(x,k,t)\ 2 , (11) 

gives the probability density of finding the particle at the point x at time t 
when initially it was on the left side of the shutter i.e., with x < and had 
a momentum k. From (HI) we see that 



\M{x^t)\ 2 = \ 



\~C( W ) 



\ ~ S[W) 



(12) 



and it is identical to the expression [4] for the intensity of light in the Fresnel 
diffraction by a straight edge. The variable w has though a very different 
meaning from the optical problem as it is now a function of time given by 
(|5[). Thus the original paper [TJ was given the name "Diffraction in time". 

All what we said above has been well known for a very long time, and 
had many applications among which we wish to mention those related to the 
time-energy uncertainty relations and decay problems 0. 

The reason that we return to this subject is that now we wish to see its 
behavior when formulated in terms of Wigner distributions functions 0] and 
also in relation with the tomographic probability developed recently by one 
of us (V.M.) and his collaborators M. 



2 Diffraction in time in Wigner distributions 
space 

Normally quantum mechanics is discussed in configuration space or, in some 
cases, in momentum space, but not in both together. Wigner 0] found that 
this limitation interfered with the application of quantum mechanics to the 
statistical physics where the description is usually given in phase space. Thus 
he introduced his concept of Wigner distributions which allow us to discuss 
some features of quantum mechanics in phase space. 

Our objective will be to formulate the diffraction in time problem, dis- 
cussed in the previous section, in terms of Wigner distribution functions. In 
this way we can visualize the phenomena in phase space and more easily de- 
termine its classical limit, and compare it with our intuitive understanding 
of the behavior of a beam of particles of definite momentum impinging on a 
shutter, when the latter is opened. 

In units in which h and the mass m of the particle are unity, and where 
the configuration space wave function is denoted by ip(x), and the momentum 
by p, the Wigner distribution function is defined as |7j 

1 r°° 
W(x,p) = (-) / ip*(x + y)ip(x - y) exp(2ipy)dy, (1) 

71 J-oo 

which has the obvious property that 

W(x,p)dp= |^(x)| 2 , (2) 



where the right hand side is the probability density at the point x, while an 
integration with respect x gives us the usual probability density [0] at the 
momentum value p. 

If we now wish to discuss the diffraction in time problem in terms of 
Wigner distributions we have to replace in ([]]) if)(x) by M(x, k,t) of (f|). 

While for our analysis k is real, we shall assume for the moment that k 
is complex with a small negative imaginary part. In this way we can use the 
expression (H) for M(x, k, t) and substituting in flU) we get 

W(x,p; k,t) = 

oo exp < — i[n(x + y) — IkH] > exp <^ i[n'(x — y) — \k'H] 



4-7T 3 J-oo J-oo J-oo (k — k*)(n' — k) 



xe i2py dKdK'dy, (3) 

where we now added the momentum k and time t to our Wigner function 
on the left hand side, as these parameters also appear in the M(x, k,t). We 
also indicate by k* the complex conjugate of k. 

The evaluation of the triple integral (§) is done in Appendix A and it 
leads to the simple result 

W(x,p;k,t) = — - — — sin S^ipt - x)(k - p^Bipt - x), (4) 

where 9 is the step function defined in (^[). Because of the presence of the sine 
function in @) we see that the Wigner distribution for the diffraction in time 
problem oscillates between positive and negative values, where the physical 
significance of the latter is not clear. On the other hand, the presence of 9 
indicates that the probability density in phase space vanishes when x > pt. 
As in our units h = m = 1, the momentum p is the same as the velocity, and 
this result is intuitively expected as the particles in our beam with momentum 
p could not yet have reached the point x. 

What is particularly interesting to us is the classical limit W(x,p; k, t) 
which is achieved when the Planck constant h — > 0. We have then to aban- 
don units in which h and m were taken as 1 and instead use standard cgs 
ones. The modifications in the form of Eq.(|]) are trivial and the resulting 
distribution function now has the form 

where 

»-!(£-)■ <«> 

If we take the limit h —> then g — > +oo as the step function takes the 
value 1 only if (pt/m) — x > 0. We can then use one of the definitions of the 
5 function 

S(h - p) = 1Sm *MLzM, (7 ) 

to write the classical limit of our distribution function as 

1 kt 

W d (x,p;k,t) = -5(k-p)9(--x), (8) 



where we used the presence of the S(k — p) in (||) to replace in the step 
function the p by k. 

We now see that our classical limit is what we expect as the only value 
possible for the momentum of the particle is p = k and besides this value is 
taken only when x < (kt/m), as for x > (kt/m) the particles would not have 
yet arrived at the point x. Thus the classical limit of the Wigner distribution 
function for our diffraction in time problem confirms our intuition. 

3 Diffraction in time in terms of the tomo- 
graphic probabilities 

In ordinary quantum mechanics the essential concept is the wave function 
which in configuration space is denoted by ip(x). From this concept one de- 
rives the probability density \ip(x) | 2 of finding the particle at point x and also, 
through appropriate transforms of ip (x) , the probabilities for given values of 
any other observables. 

Recently a change of emphasis has been proposed in which the central 
concept is the probability itself, but defining it in a tomographic way |5], [TO, 



11 . This allows us to analyze thru a single concept the probability either 



in configuration or momentum space as well as for variables that are linear 
combinations of both. The tomographic probability density |8], |TTJ was given 
in terms of the Wigner distribution through the transform 

1 /*oo /*oo /*oo 

w(X,fx,u) = — / / W(x,p)e~ lz{x - Xfl - p ' /) dzdxdp, (1) 

ZlT J— oo J— oo J— co 



where X is the position considered in an ensemble of references frames [|| [IT 
which are rotated and scaled with respect to the initial ones through the 
parameters /i, v. As an example we have that when \x = 1, v = the X 
corresponds to the normal position coordinate, but when \x = 0, v — 1, it is 
related with the momentum observable. 

In ([[]) the W(x,p) is the Wigner function defined in (|TJ) and substituting 
it in (JJ) the tomographic probability density w(X,fi,i>) is given in term of 
the configuration wave function ip(x) by 

I /*oo /*oo /*oo /*oo 

w{X,ii,v) = — / / / i){x-y)4)*{x+y)e %2py e- izi - x - x »- pv) dzdxdpdy 

ZlT J —oo J —oo J—oo J —oo 

(2) 



The integration with respect to p gives us the expression 



(3) 



and substituting it in (§), and carrying out the integration with respect to 
y, we obtain 

1 rCO fCO y-ii y-ij 

w{X,n,u) = — / ^ x + -W{x--)e-^ x -^dzdx (4) 

Z7T J -co J -co £ Z* 



Introducing now the variables 



zv 



zv 



X 



(5) 



we see the volume element dzdx in (f|) becomes drdu/\v\, so in terms of 
w(X,/i, v 

w(X, n, v 



u,r,w(X,ji,i>) becomes |IU | 
1 



2ir\u\ 



co rco 



i/j(u)ip*(r) exp 



2n\v\ 



oo J — oo 
2 



.u — r 



v 



x — /i 



r + u 



x(X,fx,u) 



drdu 
(6) 



where 



/oo 
^(u)exp[i 
-oo 



/i 2 X 

— w — u — 
2v v 



du. 



(7) 



Thus, contrary to the Wigner distribution function, the tomographic 
probability density is always positive definite. 

We now turn to the problem of diffraction in time which means replacing u 
by x in (|5]) and then ip(x) by M(x, k, t) given in terms of its expression (TO) 
containing the Green function of the free particle motion. The expression 
x{x, ft, v) takes then the form 



x(X,fi,v) 



00 f° exp[i(x - x') 2 /2t] 



oo J — oo 



e e " e ^ dxdx 



(8) 



This integral is evaluated in a straight forward but laborious way in Ap- 
pendix B, where its value is given. As we are only interested in its absolute 



value squared multiplied by (27r|z/|) l which, from ( 
graphical probability density, we see that it becomes 



gives us the tomo- 



w(X, p, v) 



2\p\ 



+ C{p) 



l 2 



+ 



1 I 2 



where 



9 



k(pt + u)-X 



(9) 



(10) 



'2n(fj,t + v) 

and C, S are the Fresnel integrals defined in Eq. (|7|) . 

We proceed now to discuss the meaning of the tomographical probability 
density given in (|8|). We mentioned above that p, v represent a rotation and 
scaling of an ensemble of reference frames in phase space with respect to the 
original one. Thus we can express them as 



/i = e r cos 9, 



v 



e T sin 9, 



(11) 



with r, 9 in the intervals — oo < r < oo, < 9 < 2tt. These expressions of 
p, v imply that our new coordinate and momenta, which we designate by 
capital X, P, are given in terms of the original ones, which we denote by 



lower case letters x,p, through the relation [12| 



x ) 


( e T cos 9 


e T sin 9 \ 


( X 


PJ- 


\ — e T sin 9 


e~ T cos 9 J 


\p 



(12) 



which is a linear canonical transformation as the determinant of the matrix is 
1. The w (X, fi, v) of @, with p given by (0); gives then the probability den- 
sity for the diffraction in time problem in the new configuration coordinate 
X defined in (0). 

If we want to return to our original configuration space we see from (|12]) 
that we must take there r = 9 = which implies X = x and /i = 1, v = 0. 



In that case p of (0) becomes 



kt 



x 



2t 



-w, 



(13) 



where w was defined in fl5|). As the Fresnel integrals are odd functions of the 
argument we have from ( |T3"D that 



C(p) = -C(w), S(p) = -S(w), 



(14) 



and thus the particular tomographic density w(x, 1, 0) becomes 

w(x, 1, 0) = \{[\- C(w)] 2 + [i - S(u>)] 2 }, (15) 



which is identical to the expression fll2|) , as we should expect. 

Thus we see that the analysis of diffraction in time phenomena in terms 
of the tomographic probabilities, allows us to study the phenomena in a wide 
ensemble of reference frames in phase space as indicated in Eqs. (9-11). 

This ensemble includes of course the original phase space (x, p) in which 
the result is given by ([TJ|) agreeing exactly with the initial analysis of the 
problem 0. 

4 Conclusion 

In the present paper the problem of diffraction in time was visualized from 
three different viewpoints. The first was the original one [|I|, in which both 
the initial conditions and the solution of the problem were analyzed in this 
same frame of reference. The solution (|5|) was given in terms of the Fresnel 
integrals, and using Cornu spiral we showed that the usual diffraction pattern 
appeared as function of time. 

In the second approach we translated our solution to the Wigner distri- 
bution space. The final expression for the probability density in phase space 
turned out to be very simple but, unfortunately, it could take both positive 
and negative values, which made its interpretation ambiguous. 

Fortunately it was possible to consider its classical limit by taking h — > 0, 
and the resulting expression (§) agreed entirely with our intuitive view, i.e., 
the probability in phase space was only different from zero when p = k and 
x < (pt/m), where all the observables are in cgs units. 

The third approach implied formulating our solution in terms of tomo- 
graphic probabilities. The latter have been introduced recently [8,10,11] to 
allow us to express the solutions in any reference frame that is rotated and 
scaled with respect to original one. In effect it implies carrying out a canoni- 
cal transformation on the original solution of the diffraction in time problem. 
Our tomographic probability solutions (^) is again expressed in terms of Fres- 
nel integrals but of an argument quite different from the one appearing (|5|). 
If the canonical transformations is the unit one i.e., X = x, P = p, then 

10 



our tomographic probability reduces to the solution fl|]), providing us with a 
check of the analysis developed in section 3. 

We finally wish to indicate that the diffraction in time phenomena de- 
rived theoretically in 1952 was, in a somewhat changed form |13|| , measured 



experimentally in 1996. Possibly a similar fate, in a distant future, awaits 
the reformulation of the phenomena presented in this paper. 
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Appendix A: Determination the Wigner function W(x,p;k,t). 
We start with the expression (M for W(x,p; k,t) and rewrite it as 



W(x,p; k,t) 



Alt 3 



OO POO 



oo J — OO 



exp[2ipy(p 



|exp[— i(nx — ^k t)} exp[i(n'x — -^k t)\ 



(«- fc*) 



(«/ - fc) 



K + ft' 

2 
dudn'. 



dy 



The first integral obviously gives the 5 function 2ir5(p 



k-\-k' ' 



(A.1) 



and so 



introducing it in (A.l) and integrating with respect to k' we obtain 



W(x,p; k,t) = 

As we have that 
1 



1 



2vr 



exp[2ip(x — pt)] 



30 exp[— 2in(x — pt)} 
oo (k- k*)(n + k-2p) 



1 



1 



1 



(k- k*)(n + k-2p) (k + k*-2p) 
we see that by introducing it in (A. 2) we get 

exp[2ip(x — pt)} 



K — k*) (« + k - 2p) 



dn. (A.2) 



(A3) 



W(x,p; h,t) 



(27r) 2 (k + k* -2p) 

J r°° exp[— 2in(x — pt)} r°° exp[— 2{k(x — pt)] 1 
'iV-oo ^ U-k*) J-oo K (n + k-2p) ]' 



(AA) 



We now note, as indicated in the text after Eq@, that we start by assuming 
that k has as a small negative imaginary part so that 



k — » k — ie; k* — » k + ie; —k + 2p — > —A; + 2p + ie 



(A5) 



Thus the singularity in the integrals in (A. 4) is in the upper half of the k 
plane. 

We can close the contour in (A. 4) by a large circle in the upper half of 
the complex k plane if x — pt < thus getting the residues of the integrals 
at the points k + ie, —k + 2p + ie. On the other hand if x — pt > we have to 
close the contour by a large circle in the lower half plane and, as the function 
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is analytic inside the contour, the integral vanishes. Then, passing to the 
limit when e — > 0, as required for our problem where k is real, we get, after 
carrying some of the multiplications, that 

W(x,p; k, t) = - — — -< exp[2i(k — p)(x —pt)\ — exp[— 2i(k —p)(x —pt)\ >. 

(A6) 
where 9 is the step function (^). As the curly bracket derived by 2z is a sine 
function we then obtain the expression (|j). 
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Appendix B: Determination of the tomographic probability w(X,fi,u) 



The tomographic probability is proportional to the absolute square of 
x{X, /i, v) where the latter is given by (0) and we rewrite it in the form 

( exp[ikx' + i(x' 2 /2t)] r 00 



where 



y/2-nit 



a 1 

a= — H . 

2v It 



/ exp[z(ax 2 — bx)]dx >dx', (B.l) 

J — oo J 



, X x' 

V t 



(B.2) 



We can rewrite the expressions in the last round bracket in (B.l) getting 

b ^2 ,ui 



ax — bx = [ y/ax — 



2^E 



6*\ 
4a)' 



As 6 2 /4a depends on x' but not on x, we first evaluate the integral 

b \2i 



exp 



y/a: 



ax 



2v^ 



dx 



11 .iff/4 



to obtain 



3 i7r/4 



z(^ + l} 



exp[i&x' + i(x' /2t)] exp 



2(^ + 1) 



e -/4 exp [^(M + l)-ll ,0 



where 



o 



i(f + 1) 



2(f + 1)' 
Using again the relation (B.3) we get 



exp[i(ax' + f3x')]dx', 



j3 = ■ " 



'fit 



+ 1) 



(B.3) 



(BA) 



dx' 



(B.5) 



(B.6) 



ax 



/2 , nl I I— I . P \2 P 



+ fix' = (y/ax' + 



2y/a J 4a' 



(5.7) 



and as /3 2 /4a is independent of x' we need to consider first the integral 



^'+5^)V 



2v^ 



;„ 2 dy 



,'.y 



,«/ 



dy 



a \/oc Jo 



2y^ 



cos|/ 2 +isiny 2 )(iy 
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v^F (1 + z) 1 v^F t 



[cf-4=)+«5(-4=) 



L 

2^' 



2^/a y/2 v 7 "^ \2^fa 

where C, S are the Fresnel integrals of (1.7) and a, (3 are given by (B.6) 
Using (B.7) to introduce (B.8) in (B.5) we obtain 



5.8) 



x(X,(i,v) 



2 ( *-■■■ 



in/A 



===== exp(-t(3 2 /4a) exp - %{X 2 I2v 2 ) 



a 1 
-v + l 



2 2^/a 



+ i 



5 + S( 2^' 



Finally replacing at,/3 by their values (B.6) we get. 



J7re in/i 
X(X, [i, v) = r —— exp 



W") 



i{X 2 /2v 2 



[i 1 



v t 



exp{—ip 



where 



P 



\ + C{p)]+i[\ + S{p) 



k{fd + v)—X 



2/i(/it + v) 



-i 



(5.9) 



(5.10) 



(5.11) 



When taking the absolute square value of x(X, /■*, v) mainly the curly 
bracket remains and thus we get Eq. (§) whose properties are discussed in the 
main text. 
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